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ABSTRACT 

A critical discussion of our knowledge about extragalactic cosmic rays and magnetic fields is attempted. What 
do we know for sure? What are our prejudices? How do we confront our models with the observations? How can we 
assess the uncertainties in our modeling and in our observations? Unfortunately, perfect answers to these questions 
can not be given. Instead, I describe efforts I am involved in to gain reliable information about relativistic particles 
and magnetic fields in extragalactic space. 
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I. What do we know? 

We know that cosmic rays and magnetic fields exist 
in clusters of galaxies for several reasons. First, in many 
galaxy clusters we observe the so-called cluster radio 
halos with a spatial distribution which is very similar to 
that of the intra-cluster gas observed in X-rays. These 
radio halos are produced by radio-synchrotron emit- 
ting relativistic electrons (cosmic ray electrons = CRe) 
spiraling in magnetic fields. We do not have direct ev- 
idence of cosmic ray protons (CRp) probably due to 
their much weaker radiative interactions. However, in 
our own Galaxy the CRp energy density outnumbers 
the CRe energy density by two orders of magnitude, 
which makes the assumption of a CRp population in 
galaxy clusters very plausible. 

Second, the Faraday rotation of linearly polarized ra- 
dio emission traversing the intra-cluster medium (ICM) 
proves independently the existence of intracluster mag- 
netic fields. It has been debated, if the magnetic fields 
seen by the Faraday effect exist on cluster scales in 
the ICM, or in a mixing layer around the radio plasma 
which emits the polarized emission (Bicknell et al. 1990, 
Rudnick & Blundell 1990). However, there is no valid 
indication of a source local Faraday effect in the dis- 
cussed cases (Enfilin et al. 2003), and the Faraday rota- 
tion signal excess of radio sources behind clusters com- 
pared to a field control sample strongly supports the 
existence of strong magnetic fields in the wider ICM 
(Clarke et al. 2001, Johnston-Hollitt & Ekers 2004). 
The detailed mapping of the Faraday effect of extended 
radio sources reveals that the ICM magnetic fields are 
turbulent, with power on a variety of scales, and with 
a power-spectrum which is Kolmogoroff-type (see Vogt 
& Enfilin, these proceedings). All these Faraday rota- 
tion measurements support magnetic field strengths of 
the order of several ^G. 

The existence of ICM magnetic fields and cosmic 
rays is not too surprising, since there are plenty of en- 
ergy sources available, which could have contributed: 



• cluster mergers: shock waves and turbulence 
{e.g., Miniati et al. 2000, 2001), 

• active galactic nuclei {e.g., Enfilin et al. 1997, 
1998), 

• supernovae remnants {e.g., Volk et al. 1996), 

• galactic wakes {e.g., Jaffe 1980, Roland 1981, 
Ruzmaikin et al. 1989), 

• decaying/annihilating dark matter particles 
{e.g., Colafrancesco & Mele 2001, Boehm et al. 
2004). 

Furthermore, we have a basic understanding of relevant 
physical processes, such as 

• particle cooling/radiation mechanisms, 

• particle acceleration by shocks and turbulence, 

• magneto-hydrodynamics. 

In principle, we just have to put all these pieces to- 
gether in order to get an understanding of magnetic 
fields and cosmic rays in galaxy clusters. This is done 
for the cosmic rays in a very schematic way in Fig. ^ 
The problem with such a generic approach is, that there 
are large uncertainties in the theoretical description of 
the different plasma processes which do not allow to 
decide a priori which are the dominant processes are 
within the complex network. Here, we will discuss two 
possibilities: 

• the relativistic electrons seen in cluster radio 
halos are re-accelerated by cluster turbulence 
{e.g., Jaffe 1980, Roland 1981, Giovannini et 
al. 1993, Brunetti et al. 2001, 2004, Gitti et al. 
2001) 
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• the radio emitting CRe are secondaries from 
hadronic interactions of a long-living CRp popu- 
lation within the ICM gas. (e.g., Dennison 1980, 
Vestrand 1982, Blasi & Colafrancesco 1999, Dolag 
& Enfilin 2000, Pfrommer & Enfilin 2004a) 

In order to have a chance to discriminate between 
the different theoretical models, detailed information 
on the electron spectrum would be very helpful, since 
the different scenarios exhibit different characteristics. 
An attempt to collect such information is displayed 
in Fig. From the figure and the discussion in the 
caption it becomes clear that the interpretation of the 
Extreme Ultraviolet (EUV) and High Energy X-ray 
(HEX) flux reported for Coma as Inverse Compton 
scattered CMB photons would require magnetic fields 
of 1.4 /xG (EUV) or significantly less (HEX). Such weak 
field strengths are an obvious contradiction to the Fara- 
day method based estimates of several /xG, especially 
if one compares magnetic energy densities or electron 
spectra as plotted in Fig. |2| 

II. What are our prejudices? 

(a) Typical cluster magnetic field strength 

Depending whom one asks, or which paper on clus- 
ter physics one consults, the assumption of the mag- 
netic field strength differ significantly. Sub-micro- 
Gauss fields are assumed when an emphasize is given to 
the Inverse Compton scattering of CMB photons into 
the EUV and HEX energy bands. Super-micro-Gauss 
fields are assumed in the context of the interpretation 
of Faraday rotation signals. 

Both methods of field estimates have their weak 
points. The inverse Compton argumentation can only 
provide strict lower limits to magnetic fields strength. 
The used observed EUV or HEX flux could have (par- 
tially) resulted from a different source or could be a 
measurement artefact. Therefore, the number of rela- 
tivistic electrons could be smaller than assumed in the 
estimate, requiring stronger magnetic fields in order to 
provide the same amount of observed synchrotron emis- 
sion. 

Faraday rotation based field estimates are also not 
straightforward, since magnetic field reversals along the 
line of sight partially cancel each other's Faraday sig- 
nal. What is left as a typical cluster Faraday signal is 
the result of a random walk in rotation measure (RM) 
in the case of turbulent magnetic fields. The statisti- 
cal RM signal depends on the statistical magnetic field 
strength times the square root of the magnetic autocor- 
relation length. The latter is unknown, and thus the 
Faraday based field estimates suffer from this uncer- 
tainty. However, the statistical properties of Faraday 
maps may allow to measure the magnetic autocorre- 
lation length under relative reasonable assumptions of 
statistical isotropy and homogeneity of the magnetic 
fields (see Vogt & Enfilin, these proceedings). 



It can be seen, that the Faraday rotation based 
field estimates can not be accommodated by sub-micro- 
Gauss field strength. For a typical cluster like Coma, 
3 — 5 ixG are reproducing the Faraday signal if a mag- 
netic length-scale of 10 kpc is assumed. Lowering the 
magnetic field strength by one order of magnitude - as 
suggested by the IC based field estimates (in the case 
of the HEX excess) - would require an increase of the 
magnetic length scale to ~ Mpc in order to reproduce 
the Faraday signal strength. But fields ordered on the 
cluster size would produce a nearly homogeneous RM 
signal, and not exhibit the many sign reversals observed 
in clusters. 

(b) The origin of radio halos 

Also the origin of the relativistic electrons produc- 
ing the cluster radio halos is controversially debated. 
There are the re-acceleration models, which are fa- 
vored by some authors since they are able to repro- 
duce nearly every observational feature reported so far. 
And then there arc the hadronic models, which are fa- 
vored by others not only because of their simplicity, but 
also since the necessary cosmic ray proton population 
should be present in the ICM because most cosmic ray 
acceleration mechanisms are believed to preferentially 
accelerate protons. 

III. How to verify models? 

It is often claimed that a model or theory is veri- 
fied. However, theories can only be falsified. To be a 
scientific theory, it has to be falsifiable (Popper). This 
means, it must be possible to derive from it unambigu- 
ous predictions for doable experiments such that, were 
contrary results be found, at least one premise of the 
theory would have proven not to apply to nature. 

Let us compare the different models for radio halos 
in the light of this statement: 

(a) Hadronic model 

The hadronic model is therefore a very good the- 
ory in the sense of Popper, since it makes a number of 
testable predictions. 

1. The hadronic model predicts gamma and neu- 
trino fluxes from clusters of galaxies. The gamma 
ray flux should be detectable with GLAST 
(Pfrommer & Enfilin and also Reimer, these pro- 
ceedings) , 

2. The hadronic model can not explain very strong 
spectral bending in the radio, since even a mono- 
energetic proton population injects a broad elec- 
tron spectrum, which produces an even broader 
synchrotron spectrum. 

3. Furthermore, the necessary energy budget of the 
hadronic model can exceed the available energy 
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Fig. 1. — Basic theory building blocks covering most of the proposed scenarios for the production and maintenance 
of relativistic particle populations in galaxy clusters and their observational consequences. Nearly all theoretical 
scenarios for cluster radio halos are summarized within this figure (e.g., re-acceleration model, hadronic model, 
hybrid hadronic-re-acceleration model, and many more). Red boxes: the main energy sources, yellow boxes: particle 
populations (CRp, CRe), green boxes: observables, white boxes: physical processes/phenomena which are (nearly) 
unobservable. The red cross indicates that the possibility that the reported non-thermal X-rays are generated by 
Bremsstrahlung of a super-thermal electron population runs into energetic problems due to strong Coulomb losses 
of such electrons (Petrosian 2001). 



sources, especially in the outer, low-density re- 
gions of clusters, where targets for an efficient 
CRp to CRe conversion are rare. 

Some of these predictions seem to put already some 
stress on the hadronic model. 

There is a strong spectral break reported for the to- 
tal flux of the Coma cluster (Thierbach et al. 2003), 
which would be too strong for the hadronic model. Ac- 
tually, this spectral break is too strong for any kind of 
synchrotron emission, indicating that there might be 
some problems with the spectral dataset of the Coma 
radio halo compiled from the literature. The strong 
spectral break might be partly explained by a con- 
tamination of the high frequency measurements by the 
Sunyaev-Zeldovich effect, however, it could also be an 
artefact of the very difficult radio astronomical mea- 



surement of extended radio halo luminosities in the 
vicinity of very strong point sources. 

It was shown that at least the radio halos in the 
Coma and Perseus clusters do not violate the energy 
constraints (Pfrommer & Enfilin 2004b). In order to 
produce the observed radio emission without too much 
energy in CRp, the hadronic models favor larger mag- 
netic field strength, in the range which are supported 
by the Faraday rotation measurements. 

The prediction of gamma rays could have not 
been tested by existing telescope sensitivities. Future 
gamma ray telescopes have the potential to ultimately 
refute the hadronic model, since there is an unavoidable 
minimum gamma ray flux of the order of the total ra- 
dio flux predicted for this class of models. If magnetic 
field strength are as low as 3 fj,G or less, the predicted 
gamma ray flux of the hadronic model has to be larger 
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Fig. 2. — Electron spectrum in the Coma cluster of galaxies (adopted from Enfilin & Biermann 1998). Red lines 
indicate parts of the electron spectrum we know to exist: thermal electrons (left), and radio-synchrotron electrons 
(right). The exact spectral location of the latter depends on the unknown ICM magnetic field strength. The upper 
limits result from physical processes, which would have produced photons in excess of observational constraints, 
if more than the indicated electrons of a given energy were present in the ICM. One process is Inverse Compton 
scattering of known background photon fields into observational bands. The other is Bremsstrahlung from energetic 
electrons with the ICM gas. Some of the used observations are actual claims of detections, and therfore could be 
measurements of the electron spectrum. Green points indicate possible explanations of the reported high energy 
X-ray (HEX) excess measured by Beppo-Sax (Fusco-Femiano et al. 1996, 2004, but see Rosetti & Molendi 2004, 
who do not find a significant signal) . Blue points indicate possible explanations of the reported extreme ultraviolet 
excess measured by EUVE {e.g., Liu et al. 1996). However, some of the marked possibilities are unlikely (marked 
by pink crosses): severe Coulomb losses would remove super-thermal and mildly relativistic electron populations 
too quickly to be replenished (Petrosian 2001, Enfilin, Lieu, Biermann 1999). The IC scattering of optical photons 
into the HEX band would require too many relativistic electrons to be consistent with other upper limits resulting 
from Bremsstrahlung. 



than this minimum prediction. 

(b) Re-acceleration models 

The re-acceleration model has not made unique pre- 
dictions which have the potential to allow to refute the 
model observationally. The model seems to be able to 
fit any yet reported radio profile and spectra due to its 
many poorly constrained parameters. Therefore, dis- 
tinctive predictions of the re-acceleration model would 
be very important, if they can be made at all (see 
Brunetti in these proceedings for first predictions of 
the re-acceleration model). 

IV. How to access uncertainties? 

Reliable spectral information on radio halos has the 
potential to refute the hadronic model. However, there 



are many known, and maybe some unknown sources 
of artefacts in radio astronomy, which can perturb the 
measured spectrum. To understand the meaning and 
significance of features and spectra we need an end-to- 
end analysis of the data reduction process. This anal- 
ysis should go from the detector signal, through (self)- 
calibration, to the map making, point source removal 
and halo flux integration step. Without such an anal- 
ysis it would be premature to claim that the hadronic 
model is ruled out due to the reported spectral steep- 
ening of the Coma radio halo. 

Also for the Faraday rotation measurements such an 
analysis would be highly desirable, in order to under- 
stand if and how observational artefacts influence our 
field estimates. In order to go into this direction, meth- 
ods to quantify the level of noise and artefacts in Fara- 
day maps were developed (Enfilin et al. 2003). They 
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Fig. 3. — Magnetic power spectrum within the cool- 
ing flow region of the Hydra A cluster of galaxies (for 
details see Vogt & Enfilin, these proceedings) . The last 
data-point at large k-valucs is likely affected by map 
noise and should be discarded. 



were used to verify the improved quality of Faraday 
maps and even the accuracy of Faraday error maps gen- 
erated with the new Polarization Angle Correcting ro- 
tation Measure Analysis (PACMAN) algorithm (Dolag 
et al. 2004, Vogt et al. 2004). These maps were then 
analyzed with a maximum likelihood power spectrum 
estimator (Vogt & Enfilin, submitted), which is based 
on the cross correlation of Faraday signals in pixel pairs, 
as expected for a given magnetic power spectrum and 
galaxy cluster geometry (Enfilin & Vogt 2003, Vogt & 
Enfilin 2003). 

The result of this exercise is not only a magnetic 
power spectrum, which is corrected for the complicated 
geometry of the used radio galaxy and of the Faraday 
screen, but also an assessment of the errors, and even 
the cross correlation of the errors. The power spectrum 
of the Hydra A cluster cool core region, which is dis- 
played in Fig. |3 exhibits a Kolmogoroff-like power law 
on small scales, a concentration of magnetic power on 
a scale of 3 kpc (the magnetic auto-correlation length) 
and a total field strength of 7 ± 2 /iG. The given error 
is the systematic error due to uncertainties in the Fara- 
day screen geometry. The statistical error is lower by 
one order of magnitude. 

V. Is a consistent picture possible? 

(a) Combining observational information 

Here, I am attempting to draw a consistent picture, 
which may explain at least a significant subset of the 

observational information* : 



• the Faraday rotation observations, which point 
towards turbulent fields strength of several /zG 
strength: in the cool core region of the Hydra A 
cluster a field strength of 7/jiG correlated on 3 kpc; 
in non-cooling flow clusters like Coma somewhat 
lower fields (say 3^iG) with a somewhat larger 
correlation length (say fO-30 kpc). 

• the radio halo synchrotron emission of CRe in 
Coma 

• the EUV excess of the Coma cluster, which may 
be understood as being inverse Compton scat- 
tered CMB light 

The Faraday measurements provide us with volume av- 

eraged^ magnetic energy densities, since the RM dis- 
persion scales as 



(RM 2 ) cx (B 2 ) V oi. 



(1) 



The synchrotron emission is also (approximately) pro- 
portional to the magnetic energy density, but weighted 
with the CRe population around 10 GeV: 



iradio OC (B 2 riCRc)vol- 



(2) 



Finally, the inverse Compton flux is a direct measure- 
ment of the number density of CRe (of the appropriate 
energy, see Fig. 



Lie oc (n C R 



c/Vol- 



(3) 



Combining the latter two measurements provides a 
magnetic field estimate (for a given or assumed elec- 
tron spectral slope), which is weighted with the CRe 
density: 



(B 2 
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(4) 



(b) Reconciling the discrepant field estimates 

The magnetic energy density derived from the com- 
bination of synchrotron and inverse Compton flux is at 
least one order of magnitude lower than the one de- 
rived from RM measurements. This discrepancy might 
be reconciled if there is a significant difference between 
volume and CRe weighted averages. This would require 

1. an inhomogeneous magnetic energy distribution, 

2. an inhomogeneous distribution of the CRe, 



tional artefacts of the difficile radio observation and data reduc- 
tion plus some impact of the thermal Sunyaev-Zeldovich decre- 
ment. 



* There are obviously some reported observations, which are not 
included in this subset: the HEX excess, since it's significance 
is debated; the reported spectral cut-off of the Coma radio halo 
emission, which may or may not be a combination of observa- 



t Actually, it is more an average with the thermal electron distri- 
bution. However, the latter varies slowly on scales over which 
the magnetic fields vary, so that it can be factorized out of the 
signal. 
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3. an anti-correlation between these two. 

The required conditions could be generated if there are 
physical mechanisms which 

a. produce inhomogeneous or intermittent magnetic 
fields 

b. anti-correlate the CRe density with respect to the 
magnetic energy density. 

A very natural physical mechanism, which could pro- 
vide this anti-correlation is synchrotron cooling in inho- 
mogeneous magnetic fields. In case of an injection rate 
of CRe which is un-correlated with the field strength, 
the equilibrium electron density is 

n CRo cx(B 2 + B2 MB )-i ; ( 5 ) 

where -Bqmb ~ 3.2^G describes the field strength 
equivalent to the CMB energy density. 

For illustration, we assume that only a small fraction 
Jb = 0.1 of the volume is significantly magnetized with 
a field strength of 10 /xG, and the rest with only 1 /iG. 
We will see later that fs =0.1 may be a plausible num- 
ber. The volume average would give {B 2 )\j^ = 3.3 [iG, 

whereas the CRe average gives (B 2 )^ c = 1.5 /iG. 
These numbers are in good agreement with the corre- 
sponding field estimates for the Coma cluster based on 
Faraday rotation and IC/synchrotron measurements, 
respectively. A larger ratio in magnetic field estimates 
could even be accomodated since the EUV emitting 
electrons are at energies below the synchrotron elec- 
trons. A spectral bump of an old accumulated electron 
population at these energies is therefore possible, and 
even plausible due to the minimum in the electron cool- 
ing rate at these energies. 

The hadronic model would provide a CRe injection 
rate which is not correlated with the magnetic field 
strength, as would be required by the above explana- 
tion of the discrepancy of magnetic field estimates by 
the two methods would work. In contrast to this, in the 
re-acceleration model one would expect a strong pos- 
itive correlation of CRe and magnetic field strength, 
since magnetic fields are essential for the CRe acceler- 
ation. 

(c) Turbulent magnetic dynamo theory 

It remains to be shown that there is also a natu- 
ral mechanism producing intermittent magnetic fields. 
The Kolmogoroff-like magnetic power spectrum in the 
cool core of the Hydra A cluster indicates that the 
magnetic fields are shaped and probably amplified 
by hydrodynamical turbulence (e.g., De Young 1992). 
Therefore, we have to look into the predictions of the 
theories of turbulent dynamo theories. 

It is generally found by a number of researchers 
that the non-helical turbulent dynamo saturates in 
a state with a characteristic magnetic field spectrum 



(e.g., Ruzmaikin et al. 1989, Sokolov et al. 1990, Sub- 
ramanian 1999 and many others). The effective mag- 
netic Reynolds number (including turbulent diffusion) 
reaches a critical value of R c w 20 ... 60. The magnetic 
fields should exhibit - more or less pronounced - the 
following properties: 

A. The average magnetic energy density eb is lower 
than the turbulent kinetic energy density £kin by 
£b ~ £kin R c 1 ■ 

B. The magnetic fluctuations are concentrated on a 
scale I, which is smaller than the hydrodynamical 
turbulence injection scale L by I « LRZ 1 ^ 2 . 

C. Correlations exist up to scale L, turn there into 
an anti-correlation, and quickly decay on larger 
scales. 

D. This may be understood by Zeldovich's flux rope 
model, in which magnetic ropes with diameter I 
are bent on scales of the order L. 

E. Within flux ropes, magnetic fields can be in 
equipartition with the average turbulent kinetic 
energy density. 

F. The magnetic drag of such ropes produces a hy- 
drodynamical viscosity on large scales, which is of 
the order of 4% of the turbulent diffusivity (Long- 
cope et al. 2003). 

(d) Confronting theory with observations 

Turbulent magnetic dynamo theory predicts inter- 
mittent magnetic fields, as favored by the proposed ex- 
planation of the discrepancy in the different magnetic 
field estimate methods. Let's see if the other predic- 
tions of the theory are in agreement with observations. 
We assume, that R c is in the range 20 to 60. 

A. The expected turbulent energy density in the 
Hydra A cluster core is of the order of 
(0.3 . . . 1) 10 -10 ergcm~ 3 , which corresponds to 
turbulent velocities of Wturb ~ (300 . . . 500) km/s. 
This is comparable to velocities of buoyant radio 
plasma bubbles (Enfilin & Heinz 2002), which are 
expected to stir up turbulence (e.g., Churazov et 
al. 2001). 

B. The expected turbulence injection scale in the 
Hydra A cluster core is of the order of (15 . . . 25) 
kpc, again consistent with the radio plasma of 
Hydra A being the source of turbulence. The 
dynamical connection of the radio source length 
scale and the magnetic turbulence scale would ex- 
plain why the Faraday map of Hydra A is conve- 
niently sized to show us the peak of the magnetic 
power spectrum. 
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C. This prediction is hard to test since fluctuations 
on scales larger than the radio source are not mea- 
surable due to the limited RM map size. How- 
ever, the downturn of the magnetic power spec- 
trum at small fc-values is at least in agreement 
with this prediction. 

D. Magnetic intermittency in form of flux ropes 
might have been detected as stripy patterns in 
the RM map of 3C465 (Eilek & Owen 2002). 

E. The fraction of strongly magnetized volume can 
become as small as /b = R^ 1 ~ 0.02. ..0.05, a 
value which is more extreme than what we as- 
sumed in our example for the Coma cluster. 

F. The expected hydrodynamical viscosity on large 
scales in the Hydra cluster is of the order of 
(1 . . . 4) 10 28 cm 2 /s. It is interesting to note, 
that a lower limit on the large scale viscosity of 
the comparable Perseus cluster cool core of 4 • 
10 27 cm 2 /s was estimated by Fabian et al. (2003). 
An upper limit on the viscosity in the (somewhat 
different) Coma cluster of ~ 3- 10 29 cm 2 /s was de- 
rived by Schucker et al. (2004). Both limits are 
consistent with our coarse estimate of the large 
scale viscosity and enclose it. 

VI. Conclusion 

To summarize, it is not clear if all observational in- 
formation can be fitted into a single consistent picture 
of extragalactic cosmic rays and magnetic fields. It 
might be that we have to revise some of our data-points. 
The choice of data-points which were used (or ignored) 
here for the construction of the presented physical pic- 
ture reflects the author's personal prejudices and may 
certainly be questioned. 

Nevertheless, it should have become clear that the 
existence of strong and possibly intermittent magnetic 
fields (several fiG) in galaxy clusters is strongly sup- 
ported by the recent detection of a Kolmogoroff-likc 
magnetic power spectra in the Hydra A cluster. Some 
of the discrepancies between Faraday-based and inverse 
Compton-based field estimates can be explained by ef- 
fects caused by magnetic intermittence, which is ex- 
pected from turbulent dynamo theory. 

Furthermore, it is argued that the hadronic gen- 
eration mechanism of the cluster radio halo emitting 
electrons is a viable model. This model is providing a 
number of stringent predictions (like minimal gamma 
ray fluxes, limits on spectral bending, maximal pos- 
sible radio luminosities), which allows detailed consis- 
tency tests with future sensitive measurements. The 
fact that this class of models seems to be more un- 
der pressure by some reported observations compared 
to alternatives as the also viable re-acceleration model 
may be caused by the larger number of predictions of 
the hadronic model in conjunction with underestimated 
systematic observational errors. 
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